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ABSTRACT 


The detection of an optical image in the presence of uniform back- 
ground light is based on a likelihood ratio formed of the numbers of 
photoelectrons emitted from small elements of a photoelectric surface onto 

which the image is focused. When diffraction is negligible and the surface 
\ r 

has unit quantum efficiency, this detectfo? is equipollent with the optimum 

detector of the image-forming light. Its performance is compared with that j 

of the threshold detector and that of a detector basing its decisions on the 

total number of photoelectrons from a finite area of the image. The 
\ 

illuminance of the image is postulated to have a Gaussian spatial distri- 
bution. All three detectors exhibit nearly the same reliability. 
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1. The Ideal Photoelectric Detector 

An optical detector is to decide whether a certain luminous object 
is present or not in an object plane at a distance R from its aperture. 

The aperture contains a lens focusing the object plane onto an image plane 
I at distance behind it. The image plane consists of a photoelectrically 
emissive surface divided like a mosaic into a large number of small areas 
-SA-^arrd^tke-^numbers- *n - -of ^photoelectrons-'emit-ted from these during an 
observation interval of duration T constitute the primary data on which is 
based the decision about the presence or absence of the object sought. In 
general, background light also enters the aperture, and an observer is to 
choose between two hypotheses: (H^) only background light is present, and 

(H^) besides the background light, light from the object is also incident 
on the aperture. 

Let the illuminance at point x of the image plane caused by back- 
ground light in the frequency interval v to v + dv be Jg(x,v)dv, and let 
the illuminance under hypothesis be 

J- (x, v)dv = J (x,v)dv + J n (x,v)dv, (1) 

where J g (x,v)dv is the illuminance produced by the light from the object 
sought (the T signal f ). Let n(v) be the quantum efficiency of the photo- 
electric surface at frequency v. We assume that the duration T of the 
observation interval is so much greater than the reciprocal bandwidth W **" 
of the object light that the numbers n^ of photoelectrons from the various 
elements of the image plane are statistically independent Poisson-distributed 



random variables with mean values 


n ik = ~ (n iiV = 


A i° 


(v) J k (x,v)d x dv/hv 


= M^Xj) 6A I » k = °» !» 


( 2 ) 


. -.the .image. , plane I, on which 

2 

d x is an element of integration, h is Planck’s constant, and 


v*> 


- T / 


n(v) J fc (x,v) dv/hv. 


(3) 


Ordinarily Mq(jc) is constant over the part of the surface I where the ob- 
ject is imaged. When the image sought is quasimonochromatic, the incident 
light will have been filtered to remove as much as possible of the background 
outside the frequency band of the image. 

The optimum strategy for deciding between the two hypotheses com- 
pares the statistic ^ 


g " 2 n i *n(n i;L /n i0 ) = n ± 


Jin 


M 0 ( 5i> 


(4) 


with a decision level g the system will choose if g > g^ and if 
g g^. This statistic differs from the logarithm of the likelihood ratio 
by a constant, 


- E Ki - n -m> = S - N . 

i 


— 8 '**£1 ** i 0 ^ w s 


(5) 
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where 


N = / M 

8 l £ 


(x)d x, M g (x) = M 1 (x) - M q (x), 


( 6 ) 


is the total average number of photoelectrons ejected by the light from 
the object. 

- .JEhe^vperf armance.of^the ^detector .can .be characterized by the 
false-alarm probability 


% Pr j 8 > g o* H oi 1 P 0^ S 0^ 


( 7 ) 


and the detection probability 


Q d " *rjg > g 0 | Hl j = 1 - P 1 (g Q ). 


( 8 ) 


Here ^(g) and P^(g) are the cumulative distributions of the statistic g 
under the two hypotheses. When hypothesis has prior probability £ and 
prior probability 1 - the average probability of error is 


(i-V 


p . ■ M 0 v 1 - v • 


( 9 ) 


and it is minimum when the decision level g^ is 


,J2] 


g Q = N s + *nU/(l - ?)]. 


( 10 ) 


Under the Neyman-Pearson criterion g^ is set to yield a pre-assigned 


-A- 



false-alarm probability Qq. 

The distribution functions Pg(g) and P^(g) can be derived from 
the moment-generating functions (m.g.f.'s) ^q(s) and h^(s), which are the 
Laplace transforms of the probability density functions (p.d.f.'s) of g 
and are given in the limit 6A^. -*■ 0 by^ 


h k (s) 



exp 


-sg 

e dP k (g) 



M-^x) 

Vi) 




k = 0, 1, 

( 11 ) 


the integral being taken over the entire image plane I. Because of (5), 
or as can be shown directly. 


h 1 (s) = expC-Nj h Q (s - 1) . * (12) 

If M (x) = 0 outside a finite area I’ of the image plane, the optimum 
s ^ 

detection statistic involves only emissions within that area, and there 
is a finite probability 



that no photoelectrons will be emitted there at all. The p.d.f.’s of g 
then have a delta-function at g = 0. 

If both M n (x) and M (x) are constant over I f , the statistic 
g[£n(M^/MQ)] has only integral values with a Poisson distribution under 
both hypotheses. When the Neyman-Pearson criterion is being used, a 
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desired false-alarm probability is set, and randomization may be necessary* 

The detection probability versus the mean total number N of photoelectrons 

s 

[3] 

ejected by the object light has been calculated and graphed elsewhere. 

In this paper we shall study the evaluation of the error probabil- 
ities and compare the performance of the optimum detector with two other, 
simpler detectors. One is the threshold detector, which utilizes a 
normalized form of the statistic 

s’ ■ n i 

The other detector simply counts all the electrons emitted from a finite 
area of the image plane. 

2. Relation to the Optimum Receiver 

The photoelectric detector just described represents a particular 
way of processing the light field at the aperture A of the optical system 
in order to choose between the two hypotheses. It should be compared with 
the optimum means of processing that field. We suppose that the object 
plane radiates incoherent, quasimonochromatic light uniformly distributed 
in a frequency band of width W about the central frequency v, W << v. 

The total radiance at point u of the object plane under hypothesis is 
B(u); under Hq it is zero. Under both hypotheses thermal background 
light of effective absolute temperature ST is also incident on the aperture. 
It is broadly distributed in direction and possesses a Planck distribution 
in frequency. 
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The field at the aperture" is considered classically as a scalar 

spatio-temporal complex Gaussian random process. In order to formulate 

the optimum decision strategy this field is broken up into spatial modes 

much as a purely temporal signal is treated by the Karhunen-Loeve expan- 
f4 51 

sion. 9 The fraction h of the object light in the g-th spatial mode 

£ 

is the p-th. eigenvalue of the integral equation 


h 

P 








(14) 


where A is the area of the aperture A and the kernel g(r) is the spatial 
coherence function of the object light at the aperture. It is the spatial 
Fourier transform of the radiance distribution B(u) of the object, 


B(u)exp(iku • r/R)d^u; (15) 

0 

where k = 2tt/X = 2ttc/v is the propagation constant, c = velocity of light, 
X = wavelength. In (14) 



® x = 3(0) = 

0 

is the total radiant power of the object. 

The coefficient of the mode function n (r) in this expansion is 
a random function of time that is in turn expanded in temporal modes that 
are eigenfunctions of the temporal coherence function of the object light. 
The combination of a spatial and a temporal mode is called a spatio-temporal 
mode. When, as assumed here, the object light has a rectangular spectral 


/ 


B(u)d^u 


(16) 
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density of width W, and the product WT of that bandwidth and the obser- 
vation time T is large, WT >> 1, there are approximately WT temporal 

modes each containing a fraction (WT) ^ of the object light; the rest 
r 4i 

are negligible. 

Under the assumption that WT >> 1, quantum detection theory shows 
that the optimum strategy forms the statistics-1 

U - V [n Jln(l + h jlr A/fW T) ] - JV , (17) 

p p s s 

E ~ 

where n^ is the total number of photons in the WT significant temporal 
components of the spatial mode p, j\^ is the mean total number of photons 
received at A from the object, and 

= [exp(hv/K t 5 r ) - 1] 1 ' (18) 

is the mean number of thermal photons per spatio-temporal mode of the 

aperture field, as given by the Planck law; K is Boltzmann’s constant. 

The number n has a Poisson distribution with mean values 
P 


i (n p lV 

IfopjHp = h p +^WT, 

E h p = 09 ) 

• £ ~ 

under the two hypotheses. Hypothesis is chosen if the statistic U ex- 
ceeds a decision level U^.. 
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When the aperture A is rectangular , x a^, with a x and a^ much 
greater than the width of the kernel 6 (r) of (14) , the eigenvalues h^ are 
approximately 


h = 6 6 

e x y 


B (p 6 


x 


Py 6 y )/B T’ 


E = 




6 = 
x 


XR/a , 

X 


6 = XR/a . 

y y 


( 20 ) 


The resolution lengths 6^ and 6^ in the object plane are much smaller than 
distances over which the radiance B(u) changes significantly. If the 
object is focused on the image plane, as in the detector of 51, the image 
is not visibly distorted by diffraction at the aperture, and if the 
imaging is otherwise perfect, the total illuminance it produces at point 
x of the image plane I is 


J (x) = CB(-Rx/R, ) 
s ~ ~ 1 


for some constant C, 

We can then write the statistic U as 


U 


■ £ 


n £n 
P 


M (x ) 
s ~p 

1 + — 7T~ 


M, 


0 


- jV 
s 


if we put 


M' = «yf'WT / 6 ' 6 ' = e/f'WTA/X 2 R 2 , 
0 x y 1 * 


m’(x) = «^s J s (x)| / J„(x)d 2 x 




-1 


( 21 ) 


( 22 ) 


(23) 

(24) 
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where 


6 ' = XR / a , 6' = AR /a (25) 

x l x y 1 y 

are the resolution lengths in the image plane. 

Comparing (22) with (4) we see that the optimum detection strategy 

handles the numbers n of photons in the spatial modes in much the same 

P 

way as the photoelectric detector handles the numbers of electrons ejected 
from the image plane. If we divided the image plane into insulated areas 
6^ x ,6^, and if the surface had unit quantum efficiency, the two detectors 
would be nearly the same. 

The moment-generating functions of the statistic U are, under the 
two hypotheses, 

\ f ( -sU, \ 

h 0 (s) = E[e |H 0 J 


= exp js^ - r yt / \TI Y) [1 - (1 + h jV l/t'WH ) _S ] j , 

' s p* £ s ) 

(26) 

h|(s) = g|e SU | « hg(s - 1). (27) 

The number of degrees of freedom in the object is conveniently defined by^ 
M - (E h p 2 ) - AA o /(XR) 2 - A o /5 x « y , (28) 


where A q is the effective area of the object, given by 
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- 1-1 


A 0 = b t 2 I / [ B (u)] 2 d 2 u 


(29) 


When M >> 1, the m.g.f. in (26) is approximately 


h 0 (s) 


exp Is ^ 

\ 8 


-A'MWT j jl 


r J\T B(u)1“ s ) , 

1 + — — (d u/A ^ , 

L e/fMWT B J ) 


B = 


VV 


(30) 


This result does not require the aperture to be rectangular, but follows 
from the two-dimensional counterpart of Theorem 1 of Chapter 5 of Grenander 
and Szego;^ f or an y function f (x) analytic at x = 0, 


y t f(h ) = M / f(B(u)/MB)d 2 u/A , M ». 1, 

£ E J 


(31) 


as can be shown by expanding f(h ) in powers of and using 


E »p" - <*V m / 


-m / „ (m) , N ,, 2 

3 (r, r)d r , 


(32) 


/ \ 

where is the ra-th iterated kernel of (14) and is proportional 


to the Fourier transform of [B(u)] m when M >> 1. The condition 

M=A/6 6 >>1 means that the image is little distorted by diffraction, 

o x y 5 

and if we assume perfect imaging we can use (21), (23), and (24) to write 
the m.g.f/s as 


h^(s) = expjs<yf 


;-/* 4 -[ i+ ^f] I 


! d2 £ > > 


M^(x) = Mq + M g (x), k = 0, 1. 


(33) 


-1L 


Comparing (33) and (11) we see that the optimum detection strategy yields 
the same performance as the optimum photoelectric detector when the 
quantum efficiency of the surface is n = 1 over the spectral band of the 
object light and diffraction by the aperture is negligible. 


3. Performance of the Ideal Photoelectric Detector 
a) Background Light Absent 

When there is no background light, Mq =0, no photoelectrons will 
be emitted under hypothesis Hq. Then under the Neyman-Pearson criterion 
a randomized strategy is required, with hypothesis chosen whenever any 
photoelectrons at all are observed; when none is emitted, is chosen 
with probability Qq and with probability 1 - Qq, where is the pre- 
assigned false-alarm probability. The probability of detection is then 

Q d = 1 - (1 - Q 0 )exp(-N g ), (34) 

where as in (6) N g is the mean number of photoelectrons ejected by light 
from the object. 

When as in a binary communication system the Bayes criterion is 
adopted, with equal error costs and with prior probabilities X, and (1 - £) 
for hypotheses and H^, respectively, hypothesis is chosen whenever 
any photoelectrons at all are observed. When none is emitted, is chosen 
provided the likelihood ratio exp(-N g ) exceeds the ratio Aq = £/( 1 - ^) ; 
otherwise H„ is chosen. The average probability of error is then 
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( 35 ) 


P = min[C»(l - C) exp (-N ) ] . 
e s 

When in particular the two hypotheses are equally likely, 

? e = \ ex P ( _N S ) » 5 = \ • 

These error probabilities are independent of the distribution of image 
illuminance. 

b) Background Light Present 

The performance of the three detectors we are concerned with, 
the ideal photoelectric detector, the threshold detector, and the simple 
counter, will be assessed under the postulate that the object creates an 
image having the Gaussian form 

M (x) = M (0) exp (-x^/2Z^) (37) 

s ~ s ~ ~ 

with width E. This might be the image of a circular nebula or, more im- 
portant, of a point source whose light has passed through a turbulent 
medium. When the turbulence can be pictured as occurring at one or a 
number of planes or 1 phase screens 1 , each introducing into the light a 
random phase shift represented by a two-dimensional Gaussian random 
process, the spatial coherence function 3 (r) at the aperture is simply 

multiplied by a factor p(r) that depends on the structure functions of 

r si 

the random phase shifts. J When these are quadratic functions of dis- 
tance, or when there are many screens, the factor p(r) has approximately 
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a Gaussian form, 


p(r) = exp(-r 2 /4L 2 ) , 


where L is a mean turbulence length that is the smaller, the more severe 
the turbulence. The image is the scaled Fourier transform of p(r) 3 (r) , 
and for a point source has the form given by (37) with the squared 
effective width 


v 2 2 _ 2 /Q 2 2 
I = X / 8tt L . 


If the observation time T is much longer than the characteristic fluctuation 

time of the random phase shifts, the distributions of the numbers of emitted 

photoelectrons will still be of the Poisson variety. We can therefore use 

2 

the m.g.f.'s given by (11), introducing the signal-to-noise ratio D and 
the total effective number of background photoelectrons through 


D 2 - H s /N 0 , 

N q = 2 it z 2 M 0 = TvfOTA/4-rrL 2 


(38) 


where rj is the average quantum efficiency over the spectral band of width 

W, T is the observation time, N is the mean total number of photoelectrons 

s 

ejected by the image, and is the mean number of background photons per 
mode, given by (18). We obtain from (11) the m.g.f.'s 
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h k (s) = exp | 


! -N„ f f[ 1 + D 2 exp (-x 2 /2E 2 )] k 


«fi- 


[1 + exp (-x^/ 2 E^)] S (d^x/2Tr 


! -N^ (1 + D 2 y) k 
0 


= exp | -w o J ^ + t 1 “ d + 0 2 y) S ]y 1 dy | , k 

(39) 


= 0 , 1 , 


after a change of integration variables and use of the circular symmetry 
of the integrand. The mean value and variance of the statistic g under 
each hypothesis can then be obtained by expanding the m.g.f. f s about s = 0, 


EtglHjJ = N 0 / (1 + k D 2 y)£n(l + D 2 y)y 1 dy, 


Var^liy = N Q I (1 + k D 2 y) £n 2 (1 + D 2 y)y 1 dy, k = 0, 1 . 

(40) 


It is not generally possible to calculate the p.d.f.'s of the 
detection statistic g in closed form from m.g.f.’s such as those in (11) 
and (39). The logarithmic likelihood ratio g has what is known as an 
infinitely divisible distribution, but although there is a vast mathe- 
matical literature about such distributions, mostly concerned with limit 
theorems , ^ we have been unable to find there any specific advice 
about calculating the probability density functions. In a previous 
paper a Gaussian approximation was used; it is valid only when 
Nq >> 1 and << 1. Farrell recommended approximating the proba- 
bility density function (p.d.f.) of the statistic g by a gamma distribution 
whose mean and variance match those of the true distribution as given by 
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( 40 ). 


As the m.g.f.-is the Laplace transform of the p.d.f., it is natural 

to try evaluating the inverse Laplace transforms of h^Cs) and h^(s) by 

the saddlepoint method* Daniels applied this method to finding the p.d.f. 

[121 

of the sample-mean statistic. We have used it to obtain the cumulative 

distributions directly by evaluating the contour integral 


i - p k (g> - 


/ 


o + i°° 


gs 


[1 - h-^(s)]e ds/27ris 


a - i°° 


a > 0, k « 0, 1, 


(41) 


which can also be written 


1 - P k (g> - - 


/ 


a + i°° 


gs 


h^(s)e ds/2Tris, 


a - 


a < 0, k « 0, 1, 


(42) 


by displacing the contour to the left of the origin. By an elaboration 
of the saddlepoint method. Rice produced an asymptotic expansion that is 
most suitable for numerical evaluation. 

The saddlepoint occurs at the values of s at which the phase 

4> k (s) = &n h fc (s) + gs ' (43) 

of the integrand of (42) is stationary; it is the solution of the equation 

<(s) = h^(s)/h k (s) + g = 0, (44) 
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which from (39) yields 


g = N Q J (1 + D 2 y) k S £n (1 + D 2 y)dy/y, 

0 

s = s^, k = 0, 1, (45) 

where the primes denote differention with respect to s. For the ideal 
photoelectric detector, by virtue of (12), 

= Sq 4* 1. (46) 

When g > E(g|H k ), < 0 and (42) is used; when g < E (g | H^) , > 0 and 

(41) is used. Rice’s asymptotic expansion requires the values of the deri- 
vatives of the phase at the saddlepoint, 

d m <f>/ds m = N Q 

0 

s = s^, m ^ 2, k = 0, 1; (47) 

these were evaluated by numerical integration. 

In Fig. 1 we have plotted the cumulative distributions of the 
statistic g under hypotheses and for a typical set of values of N g 
and Nq, as calculated by the saddlepoint method and from the Gaussian 
and gamma-function approximations. The latter two approximations are 
least accurate in the tails of the distributions, and for small false-alarm 
probability there may be a serious error in the decision level g^ if 


J (1 + D 2 y) k S £n m (l + D 2 y)dy/y, 
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either is used. 


Under the Neyman-Pearson criterion the false-alam probability 

Qq is pre-assigned, and it is necessary to determine the decision level 

and the associated saddlepoint Sq. An iterative search method was 

used that took advantage of the fact that as s^ becomes more negative, g^ 

increases and decreases. The behavior of the saddlepoint s^ as a 

2 

function of the signal-to-noise ratio D is depicted in Fig. 2 for a 

typical case. The probability of detection is plotted as the solid 

curves in Figs. 3 and 4 as a function of the mean number N g of received 

photons for various values of the mean number of background photons 
2 

in an area 2rr E concentric with the Gaussian image. In Fig. 5 is 
plotted the average probability P^ of error as calculated from (9) for 
the ideal detector when used in a binary communication system; £ = 

and g 0 =N s- 


4. The Threshold Detector 

The structure of the optimum detector depends on the signal-to- 
2 2 

noise ratio D , and if D is unknown, the detector must be designed for 

2 

some reasonable standard value. When D is very small, the optimum de- 
tector is nearly equipollent to the threshold detector, which is 

based on the term of lowest order in an expansion of the logarithm of 

2 

the likelihood ratio in powers of D . For detecting an image on a 
photoelectric surface, the equivalent threshold statistic is^^ 


g 0 = n i m(x) = M g (x)/D 2 M 0 (x) , 


(48) 
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which is compared with a decision level Sqq> hypothesis being chosen 

when g > g Its probability distributions under the two hypotheses 

0 0U 

can be determined from the m.g.f^s 


h ^ 9 ^(s) = E[exp(-g Q s) | h^] 


= expj- f M q(x ) [1 + D 2 m(x)] k |l - exp [-sm(x) ] j d 2 x j , 


k = 0, 1 , 


(49) 


which for our Gaussian image in (37) becomes 


h^ (s) = exp 


j-N Q J (1 + D 2 y) k (l - e Sy )y 1 dy 
0 


= exp | -N n [C + £n s - Ei(-s) + kD^(s -l + e S )s^]J, 


k = 0, 1, 


(50) 


where Ei(x) is the exponential integral and C = 0.577215... is Euler's 
constant. The derivatives needed for the asymptotic expansion resulting 
from the saddlepoint method can be written down in closed form, 

d m Un h^ e \s)]/ds m = (-l) m N qJ y™ 1 e Sy (l + D 2 y) k dy 

0 

= (-1)® N q | s~ m [ (m-1) ! -r (m, s) ] + kD 2 s _(m+1) [m!-T (m+l,s) ] 

(51) 
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in terms of the incomplete gamma-function 


/ CO 

t m e *■ dt = (m-1) ! e X x r /r!. (52) 

r = 0 
x 

The decision level g which is independent of the signal-to- 

0 0 

noise ratio, can be determined as described in Section 3, after which 
the detection probabilities are calculated by the asymptotic .series derived 
from (41) and (42), into which h| ; (s) from (50) is substituted. The 
results are plotted as dashed curves on Figs. 3 and 4 for = 0.5 and 
Nq = 5. They approach the detection-probability curves for the ideal 
photoelectric detector in an osculatory manner. 

If the image detector is to be used in a binary communication 
system, the decision level minimizing the average probability of error , 
must be found from the likelihood ratio 

dp i^ 5 * * 8 e^ dp o^ 8 e* ~ & 9 8 eo > 

and this requires finding the p.d.f.’s of g by inverse Laplace trans- 
formation of (50), after which the minimum value of P^ is determined 
from (9). We have not carried through this calculation. 


5. The Simple Counter 

The simplest way to detect the image is to count all the photo- 

electrons emitted from a finite region of the image plane that is concentric 
with the expected image illuminance J (x) . The total number n of 

g r^f 
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photoelectrons is then the datum on which the choice between hypotheses 
Hq and is based. It will have a Poisson distribution with different 
mean values under the two hypotheses, and in order to attain a pre- 
assigned false-alarm probability, randomization will in general be 
[16] 

necessary. 

Applying this scheme to detecting the Gaussian image given by (37), 
we assume that we count all the photoelectrons emitted from a circle of 
radius p£ centered at the origin. Under hypothesis the mean number 
of counts is 


E(n|H Q ) = n Q 


ttM,. 


2 

P 



1 2 „ 

2 P N 0> 


(53) 


where is given by (38). Under the mean number is 


E(n| H 


i> = n i = JJ Vs )d2 s = N o|} p2 + 1)2 < 1 - e " P /2 > ] ’ 


x <pZ 


(54) 


where D = N s /Nq is the signal-to-noise ratio. The detector chooses 
hypothesis whenever the number n of counts exceeds an integer 0; 
is chosen when n < 0; but when n = 0, is chosen with probability f. 
The false-alarm and detection probabilities are then 


Q q = exp(-n 0 ) |fn 0 6 /e! •+ £ \ / n! ’ 

n - 6 + 1 J 

Q d = exp(-n 1 ) j^f n;L e /© ! + ^ n^/nll . (55) 

n = 0 + 1 J 
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In Figs. 3 and 4 the probability of detection so calculated is plotted 

versus the mean number N g of photoelectrons ejected by the Gaussian 

—3 —5 

image, for Qq = 10 and 10 . For each point on the curve the value 

of p was varied to yield maximum detection probability. The maxima are 
quite flat, and the radius pE of the counting area is not critical. The 
values of p used are listed in Table I. 

When the average error probability in (9) is to be minimum, the 
likelihood ratio A(n) formed from the number n must be compared with the 
decision level A Q = £/(l - 4), and hypothesis H 1 is chosen when A(n) > A Q . 
Thus the image sought is declared present whenever 

n > j £n[ £/ (1 - O] + n ! “ n 0 | / £n ( n 1 /n 0 ) ; 

otherwise hypothesis is chosen. The average probability of error as 
given by (9) can then be calculated by summing Poisson probabilities, and 
it has been plotted as the dashed curves in Fig. 5 at C = ^ For each 
point on the curve the value of p was chosen to yield minimum error 
probability. 


6. Conclusions 

From Figs. 3, 4, and 5 we see that the ideal photoelectric detector 

and the threshold detector, which utilize information about the shape of 

the image J (x,v), are hardly better than the simple counter when the 
s 

effective area of this is selected to yield maximum detection probability 
or minimum error probability P^. If the image has a uniform distribution 
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over a finite area I* of the image plane, outside of which it vanishes, 
all three detectors are basically identical. Because the ideal photo- 
electric detector is nearly equivalent to the optimum detector of an 
incoherently radiating object in the presence of thermal background 
light when it has unit quantum efficiency and diffraction is insignificant, 
there is little one can do to improve the simple counter except to 
increase its quantum efficiency; registering the points of emission of 
the photoelectrons does not help much to increase the reliability of 
detecting an image. 
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TABLE I 


Optimum Radius for Simple Counter 


N 


0 


0.5 


N 

s 

p(Q 0 = 10 3 ) 

P(Q 0 = 10 5 ) 

1 

1.31 

1.37 

2 

1.72 

1.68 

4 

1.72 

1.68 

8 

2.10 

1.68 

12 

2.10 

1.98 

20 

2.10 

1.98 

24 

2.46 

2.28 


N o - 5 

N 

s 

p(Q 0 = 10" 3 ) 

p(Q 0 = 10 5 ) 

1 

1.36 

1.37 

2 

1.44 

1.44 

4 

1.60 

1.58 

8 

1.60 

1.58 

12 

1.75 

1.65 

20 

1.82 

1.78 

24 

1.82 

1.78 






FIGURE CAPTIONS 


Fig. 


Fig. 2 


Fig. 3 


Fig. A 


Fig. 5 


False-alarm probability Qq and detection probability as 
a function of decision level g^ for ideal photoelectric 
detector as calculated by the saddlepoint, Gaussian, and 
gamma approximations; Nq = N g = 5. 

. Progress of the saddlepoint Sq of the phase of the inte- 
grand of (42) as a function of the signal-to-noise ratio 
D 2 ; N = 0.5 and 5, Q q = 10 -3 and lO -5 . 

. Probability of detection of a Gaussian image vs. the 

mean number of photoelectrons ejected by the image; 

_3 

Qq = 10 . Curves are indexed by the’ noise parameter 

Nq defined in (38). 

. Probability of detection of a Gaussian image vs. the 

mean number N of photoelectrons ejected by the image; 

s 

Qq = 10 Curves are indexed by the noise parameter 
Nq defined in (38). 

. Average error probability P^ in deciding between presence 
and absence of Gaussian image when these have equal prior 
probabilities, vs. mean number N g of photoelectrons 
ejected by the image. Curves are indexed by the noise 
parameter Nq defined in (38). 
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Probability of detection of a 
Gaussian image vs. the mean number 

N of photoelectrons ejected by 

s -5 

the image; Q q = 10 . Curves are 

indexed by the noise parameter 

defined in (38) . 
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Fig. Probability of detection of a Gaussian 

image vs. the mean number N of photo- 

S -3 

electrons ejected by the image; = 10 
Curves are indexed by the noise parameter 
Nq defined in (38) . 
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Fig.^. Average error probability P e in deciding 
between presence and absence of 
Gaussian image when these have equal 
prior probabilities, vs. mean number 
N g of photoelectrons ejected by the 
image. Curves are indexed by the noise 
\ parameter N n defined in (38). 
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